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Abstract 

We give an exhaustive description of bifurcations and of the number of 
solutions of the vacuum Lichnerowicz equation with positive cosmological 
constant on x 5^ with U{1) x S'0(3)-invariant seed data. The resulting 
CMC slicings of Schwarzschild-de Sitter and Nariai are described. 
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1 Introduction 


One of the challenges of mathematical general relativity is to provide an ex¬ 
haustive description of all physically significant solutions of the constraint 
equations. In mathematical cosmology the relevant model is usually taken 
to be that of spatially compact solutions. For definiteness vacuum models 
only will be considered here. An exhaustive description of such models 
with constant mean curvature (CMC) r has been given by Isenberg [14] 
using the conformal method, assuming that the cosmological constant A 
vanishes. The analysis there carries word for word to the case 

2n 

_ 2 ) ^ ’ , ( 1 . 1 ) 

in space dimension n, where the question whether is identically vanish¬ 
ing or not is replaced by the corresponding question for It 

thus remains to understand what can be said about the conformal method 
for constructing solutions of the CMC vacuum constraint equations when 

(1.1) does not hold. 

When the extrinsic curvature tensor K is pure trace and when (I.l) 
fails, the Lichnerowicz equation for the conformal factor reduces to the 
Yamabe equation with positive scalar curvature. Already one of the sim¬ 
plest models, namely x 3“^ with the standard product metric, provides 
an example where solutions of this Yamabe problem are not unique [29]. 
It turns out that this case can be described in an exhaustive way (see [30] 
for a comprehensive and clear analysis). The object of this note is to ex¬ 
amine this same model from the point of view of the conformal method 
for constructing solutions of the general relativistic constraint equation. 

As such, we consider U{1) x S'0(3)-symmetric seed fields on x 
for the conformal method. Given a cosmological constant A which we 
assume to be positive throughout, such fields can be parameterized by the 
length T of the A^-factor in the metric, the curvature scalar .R > 0 of the 
S^-factor of the metric, the trace r of the extrinsic curvature tensor, and 
the norm of a seed TT-tensor which we encode in a parameter a; see 

(2.2) -(2.3) below: 

Theorem I.l Let R, A gM.* and a, r € M. Let (f) denote a solution of the 
Lichnerowicz equation 



on X 3“^ endowed with the product metric g defined in (2.2) and with 
seed fields as above. Assuming further that < 3A and a := L ^ 0, the 
following holds: 

1. The equation has no solutions when 



2. When the inequality in (1.2) is replaced by an equality there exists 
precisely one solution, which is constant. 
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3. When the inequality symbol > in (1-2) is changed to <, solutions 
exist and all are SO{i)-invariant. Moreover: 

(a) There exists a function T{a,T,R) > 0 (see (5.2) together with 

(4.11) and (4.15)J, with T ^ oo as ~ approaches 

from below, such that for 

T € {nT{a, r, R) ,{n l)T{a, r, R)] 

there exist exactly two constant solutions and exactly n non¬ 
constant solutions, counted modulo isometry. 

(b) Furthermore, defining as the largest integer k such that 

there exist explicit constants 

/ o \ 3/2 

1 I r\ 

I 3" ) > ai > • • • > afc„ax > 0 

such that, for each a in the range {—ak,—oik+i] U [ak+iiOik) 
(resp. (— the Lichnerowicz equation has 2 con¬ 
stant solutions and k (resp. kmax) non-constant solutions, counted 
modulo isometry. 

We give an explicit expression for the period function T of point 3(a) in 
Equation (5.2) below (compare (4.11) and (4.15)), assuming the conformal 
gauge R = where —^ := ^ — that can be achieved by rescaling 

the metric g. A general explicit expression can be obtained by solving 
(6.4) with k = \ for T as a function of a, but is not very enlightening. 

It should be clear that the case of point 2 of Theorem 1.1 is dramati¬ 
cally unstable: small perturbations of the parameters might lead to non¬ 
existence, while it follows from the analysis below that there exist small 
perturbations of the initial data for the solution that lead to a conformal 
factor which reaches zero in finite time, so that the conformally rescaled 
metric does not correspond to a complete periodic geometry. It should 
also be clear from what is said below that the smaller constant solutions 
of point 3 of Theorem 1.1 are again unstable in the sense just given, and 
that the larger constant solutions, as well as all non-constant solutions are 
stable. 

Our analysis is restricted to three space-dimensions, but we expect very 
similar results to be true on x A" with U{1) x S'0(n-|- l)-invariant seed 
fields. We note interesting dimension-dependent phenomena arising in the 
Yamabe problem [4,16], with similar behaviour expected to occur for the 
higher-dimensional Lichnerowicz equation. 

Recall that large families of non-trivial solutions of the Yamabe prob¬ 
lem can be constructed using bifurcation-theory methods, see [11,21] and 
references therein. In Section 6 we apply these methods to our case. It 
should, however, be recognized that in view of the symmetry result of [15], 
a direct analysis of the PDE applies and provides immediately a much 
clearer picture. 
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It follows from the (generalized) Birkhoff Theorem that the initial data 
resulting from our solutions of the Lichnerowicz equation can be realized 
by CMC sections of Schwarzschild-de Sitter space-time, or Nariai space- 
time. We discuss this in Section 7; compare [2,10]. 

It should be said that non-existence, or existence of multiple solu¬ 
tions, for the Lichnerowicz equation have been already pointed out in 
the literature in settings much more general than ours [1,9,12,13,17- 
19, 22, 24,31,32]. The interest of our model resides in its simplicity, which 
allows us to obtain a complete description of the set of its vacuum solutions 
using elementary arguments, drawing upon the deep results of [5,15]. 


2 The model 

The initial data manifold M we consider is x S"^. The initial data 
metrics 

9 = <l>^g (2-1) 

will be conformal to 

9 = 9f,R-= + (2.2) 

where tp is a. 27r-periodic coordinate on S'^, T and R are positive constants, 
while dfl^ is the unit round metric on S'^. The metric 9f has scalar 
curvature R. A constant rescaling of gf can be absorbed in a redefinition 
of T and R which leaves invariant the product T^R. We will write g instead 
of gf f when the explicit values of T and R are not essential. 

Following [6], the extrinsic curvature tensor K will be taken of the form 

K = j + ^-g =: + ^-g , (2.3) 


where a and r are non-negative constants. This is the general form for a 
symmetric 17(1) x S'0(3)-invariant 2-tensor. Note that the “seed tensor 
field” L is ^-transverse and traceless. The multiplicative normalisation 
factor in L has been chosen so that \L\^ = jaj. 

In the current situation the general relativistic constraint equations 
will be satisfied by {g,K) if and only if 


^g4>- 




(2.4) 


where 

^ _ A 
■“ 12 ~ 4 ’ 


(2.5) 


with /? > 0. We have introduced the notation “/3^” to emphasise the fact 
that we focus on the case where the right-hand side of (2.5) is negative. 
In fact, when 0^ is negative and 0 0 the solutions are unique (cf., 
e.g., [14]). This implies that (p inherits then the symmetries of the metric 
g, and hence is constant. We will see that this is not always the case 
anymore when is positive. 
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From now on we assume > 0 (otherwise this is the Yamabe problem, 
already discussed in the references pointed out in the introduction) and 
that > 0 (otherwise there is only a constant solution). 

It is easy to see that a positive solution of (2.4) exists if and only if a 
constant solution exists. For this note that, since Ap integrates to zero, 
there exists a point p on x 5^ such that A(j>{p) = 0. At this point we 
have 

- Y</'(P)■^ (2.6) 

which implies that the constant function (p = (j){p) solves (2.4). 

A first corollary of this, keeping in mind (2.6), is that i? < 0 is in¬ 
compatible with the existence of positive solutions of (2.4) with /3^ > 0, 
therefore 3“^ cannot be replaced by another two-dimensional surface in our 
model. 

We continue with the following result, where we allow R, a and /? not 
to be constant: 

Proposition 2.1 Consider (2.4) on a compact Riemannian manifold with 
continuous functions R, a and P satisfying P^ > 0. If 

min min/?^ > ^ max , (2-7) 

then (2-4) has no positive solutions unless R, and P'^ are all positive 
constants and the inequality in (2.7) is an equality, in which case there is 
a unique positive solution, which is constant: 



Proof: Write (2.4) as A(j> = F. A simple analysis of the polynomial 
(p — Rp^ + gives F > 0 when (2.7) holds. Multiplying the 

equation by p and integrating by parts gives p = po = const, F{Pq) = 0, 
and the result readily follows. □ 

As is well known, the problem at hand is conformally covariant. It 
follows that we can always rescale g so that a constant solution, whenever 
one exists, equals one. After such a rescaling we will obtain 

R = P‘^ + a‘^ . ( 2 . 8 ) 


This normalisation will be often used in what follows. 


3 Solutions depend only upon 'ip 

Uniqueness and non-uniqueness results for solutions of several classes of 
semilinear equations on S^, possibly with isolated singularities on the 
north and south pole, have been proved in [15]. Here we check that The¬ 
orem 1.5 there applies and shows that: 

Theorem 3.1 Solutions of our model depend at most upon p. 
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For the sake of completeness, we recall the result from [15] we will need. 
The conformal Laplacian .jSfgs on the unit three-dimensional sphere reads 

■^53 = As3 - - . (3.1) 

We denote by n = (0, 0, 0,1) G S'^ C ffi"* the north pole, by s = (0,0,0, —1) 
the south pole, and we set 

:= \ {n, s} . 

We denote by u the latitude on S^, ranging from 0 at the north pole 
to TT at the south pole. We will need the following special case of [15, 
Theorem 1.5]: 

Theorem 3.2 Let f G C°((0,7r) x (0,oo),M), and let v G (7^(12,M) be a 
positive solution to 

-^S3V = f{u,v). (3.2) 

Assume that f satisfies the following conditions: 

1. for each 9 G ft the function s i—)■ s~^f{9,s) is non-increasing on 
(0,oo), 

2. for each s G (0, oo), the function 9 i-G- f(9, s) is decreasing strictly on 
(0,7r/2) and increasing strictly on {Trl2,Tr). 

If 

liminfu(6*) > 0 and liminfu(6*) > 0, (3.3) 

9^n 0—>s 

then V is rotationally symmetric about the line through n and s (or equiv¬ 
alently, V depends only on u). 

Note that in (3.3) liminf is allowed to be infinite. 

We are ready to pass to the proof of Theorem 3.1: 

Proof of Theorem 3.1: The idea of the proof is to view (2.4) as an 
equation over M x 5”^, the universal cover of 5”^ x and map it, via a 
conformal isometry, to an equation on \ {n, s}. The round metric g 
on can be written as 


g = d9^ -\- sin^ 9dLi^ . 

We introduce the coordinate t = so that the metric g reads 

2 

g = dt^ -\- -^dLl^, 

R 

and seek functions t = t{9) and pL = p,{9) > 0 so that 
g = pf {d9'^ -\- sin^ 9dLl^) . 

Identifying the coefficients of d9 and we are led to the following 
relations: 

/ 2 \ 

dt = pfd9 and =P^An9. 

As a consequence, the functions t(9) and p are given by 
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The conformal Laplacians of g and g are related by the following well- 
known formula: 

Hence, viewing (2.4) as an equation on H and setting v = g,(j), we have 


—^gV = 



-v^ + 


,,-7 


m sin® e 


(3.4) 


The assumptions of Theorem 3.2 concerning / are readily checked. Since 
(j) is a. positive function on x S'^, it is uniformly bounded from below. 
So, g{9) —>• oo as 0 —> n or 0 —> s, i.e. 


lim inf v{9) = lim inf v{9) = oo. 

0 —^ri 0 — 

We conclude that Theorem 3.2 applies to Equation (3.4): the function v 
depends only on 9. Since /r depends also only on 9, /pis a function of 9 or, 
equivalently, a function of t. □ 


4 ODE analysis 

In view of Theorem 3.1, we seek solutions </> of the Lichnerowicz equation 
(2.4) with R given by (2.8) such that 

(j) = (j){'iP), d^(j){0)=Q-, (4.1) 

note that the last condition can always be fulfilled by an adequate choice 
of the origin on the circle. Thus 

(27r)^ _ 1 ^ / 2 I I a2j,5 I 2J,-7^ 

= + (4.2) 

The conserved energy for (4.2) reads 

^) = ^02 _ + g <^® - =: 1^2 ^ ^ (4 3 ) 

where a dot denotes a derivatives with respect to 

(4.4) 

Keeping in mind our assumptions </> > 0, 7 ^ 0 and /3 > 0, the equation 

dV/d(j) = 0 can be written as 

{y — l){x^ + x^y — 2y'^) = 0 , where y = (j)^, (4.5) 

and where we set 

a: := |a|v^//3 > 0. (4.6) 

The positive solutions are y = 1 and y = j{x + \Jx"^ + 8 ), distinct unless 
a; = 1. Representative plots of V can be found in Figure 4.1. 

When |a| = ^ j the only solution which remains bounded away from 
zero for all times is ^ = 1. This case corresponds precisely to that already 
covered in Proposition 2.1, and thus from now on we assume that 

|a| 7 ^ /S/v^ or, equivalently, a; 7 ^ 1 . 
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Figure 4.1: Typical form of the potential V{(j)) with |q;| < (left), 

|a| = /3/-\/2 (middle) and |a| > /3/-v/2 (right). 

4.1 Solutions on M x 

Let us relax the condition that ip is a periodic coordinate, and consider 
instead (4.2), where the '0-derivatives are replaced by derivatives with 
respect to the parameter t G K of (4.4). The nature of the solutions 
M 9 t I—>■ p(t) is apparent from the phase portrait in the (0,0) plane of 
Figure 4.2. 




Figure 4.2: Phase portraits with /3 = 1, a = .2 (left) and /3 = 1, a = 1.2 (right) 
in the (0, 0) plane. The color encodes the value of the energy 11(0, 0). The red 
curves correspond to the level sets of the instable constant solutions. 

Let us discuss some overall features of the solutions. 


4.2 |a|<0/\/2 

The critical point (0 = 1,0 = 0) is stable if and only if |a| < /3/-\/2, as 
should be clear from Figure 4.1. The associated critical energy is 

lli=ni) = -^(2a"+/?") • (4.7) 


The second critical point (0 = 02 < 1, 0 = 0) has energy which we will 
denote by H 2 = H2{a,j3). An analytic expression for H 2 can be obtained 
but is not very enlightening: 


H2 


x'^ -f -1-16-1- (x^ -I- 2x)-\/x^~+~8 
372 


48 


(x -I- Vx'^ + 8) 


(4.8) 


























There is an orbit corresponding to a non-trivial homoclinic solution 
with energy <j)) = H 2 which asymptotes to (j )2 as t tends both to plus 
and minus infinity. On both pictures of Figure 4.2, this solution lies on 
the piece of the red curve that closes up. 

All orbits lying in the conditionally compact set, say If C {((/>, G 
enclosed by this homoclinic orbit are periodic. These are the only orbits 
with (j) bounded and bounded away from zero, and hence the only ones of 
interest to us as solutions of the Lichnerowicz equation on x 5^ leading 
to a spatially compact vacuum data set with the same topology. 

The periodic orbits oscillate between i^min (a, /3, E) and 4>-max (a, /3, E), 
where 

denotes the energy of the solution. It should be clear from Figure 4.2 that 
the function E 1 —>■ ^min (a, /3, E) is monotonously decreasing to 0„iin («> /3) = 
(j) 2 , while E I—>■ ^max (a, /3, E) is monotonously increasing to a value 4>max (a, /3). 
There is a bound 

^max {a, P) <V3 

which is approached as a —>■ 0. It is attained on the solution with a = 0 
with energy E = H 2 , for which (/)2 = 0; this solution closes-off M x to 
a smooth round S^. A plot of {ct, P) can be found in Figure 4.3. We 
note limQf,^|y^^_^Q 0. 



Figure 4.3: («, P) and </>max («) P) as functions of the scaled variable x 

V2\a\/p. 


The discussion so far applies to solutions on K x 5"^. As such, given a 
metric on x we wish to find all solutions of the Lichnerowicz equation 
in our context. Now, a solution on M x with minimal period T leads 
to a solution of the Lichnerowicz equation on x 5^ with metric g^j, 
for any integer a > 1, by replacing the 5'^-factor by its n-fold cover. The 
question that arises is then which values of T are realised by the solutions 
above. To answer this we need to understand the period function. 


4.2.1 The period function 


Consider the function which to a periodic orbit with energy E associates 
its minimal period T{a, /3, E). For any such orbit with p varying between 
pminic^,P,E) and prasLx{.Oi, p,E) the period equals 


T = V2 




dp 


y/E-V{p) ’ 


(4.9) 
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where the turning points l3, E) and j3,E) are found by 

solving the equations 


P,E)) = E = y P, E )), 


with (j)n,in{a,/3,E)) G [(j) 2 {a, P) ,1] and E)) G [l,oo). Since in 

our case y is a real analytic function of </>, the real analytic version of 
the implicit function theorem shows that away from the critical level sets 
of E[ the functions E i—>■ (pmin and E i—>■ (/)max are real analytic; compare 
Lemma 6.6 below. 

When E approaches the energy of the stable critical point, (/)s, where 
V"{(j)s) > 0, the period approaches that of linearized oscillations around 


T 


27r 


^/V^s 


(4.10) 


In particular, when |a| < fily/2 the stable critical point is (()i = 1 and one 
has 

T^T,{a,P) = -^==. (4.11) 

Near to and away from the critical point (j) = 1 the function T is differ¬ 
entiable, with the sign of the derivative of T with respect to E determined 
by the sign of the Chicone test function [5] 


N = {G'f 



n 


(4.12) 


where G(<()) = V{(j)) — y(l) is the potential normalised so that G(l) = 0, 
on the interval P)), P-mmioc, P))]- N can be computed and takes 

the following form: 


N = 


P^iP^ - 1 )" 

768(()22 


P{P), 


where P is a polynomial of degree 28 in <j) which is conveniently computed 
with e.g. Mathematica. Setting 


p2 ■— P^ , 

one finds by inspection that the polynomial X i—)■ P{p 2 + X) in X has 
all its coefficients positive for all x G [0,1]. This is, in fact, obvious for 
all coefficients except possibly for the term linear in X and the constant 
term. Now, the coefficient of X in P{p 2 + X) equals 

——^^^(1 - (1 - p 2 )p 2 ){^ + P 2 ) {p 2 + 2 ) 

y.{l + p2 + Pi) (62 -h ASPI - 430^ -b 2401 - 670^) , (4.13) 

where we expressed x (defined in (4.6)) in terms of 02 using (4.5): 


Note that 



P2 e (0,1) 
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since x G (0,1), so 1 — (1 — (t> 2 ) 4’2 > | and we are left with proving that 
62 + 48(/)^ - + 240® - 670^ > 0. 

This follows from the following observation: 

62 + 480^ - 4301 + 240^ - 670| > (62 + 48 - 43 - 67)0^ + 240| = 240| > 0. 

Finally, the coefficient of X® can be written in the form 

480f (1 - 0®) (0® + 04 + 202 + 2)% (4.14) 

which is again manifestly positive since 02 < 1. 

From the fact that all coefficients of P(X+02) are positive, one immedi¬ 
ately concludes that N is non-negative on the interval [0min(Q;, 0)), 0max(Q:, 0))], 
thus proving that the period function is increasing with E. 

When moving continuously amongst the solutions so that their energy 
E tends to 772, the period of the solutions grows to infinity as is to be 
expected since the (bounded) solution with E = 772 is a homoclinic orbit. 

A plot of the period function E ^ T{a,l3,E) can be found in Figure 
4.4 for a = 0.2 and 0=1. In this case one has 77i ~ —0.045 and 
772 -0.0364. 



Figure 4.4: Values of the period of oscillation with a = 0.2 and 0 = 1. 


4.3 |a|>0/\/2 

The analysis is very similar to that of the case |a| < 0/-\/2- In this case, 
we have x G (l,oo). The stable point becomes 02 and a calculation shows 
that ^ 

F "( 02 ) = ^\/% 2~+~8 ^ 3 a ; - \/ a ;2 + 8 ^ 

which is clearly positive if and only if x G (l,oo). When the energy 
of a periodic solution approaches 772, its period approaches that of the 
solutions of the linearized problem around 02 : 


T^T2{a,l3) 


2tt 


4y27r 1 1 

101 (x2 + 8)'/%3x-V^^)'^’ 


(4.15) 

As expected, the period of small oscillations goes to infinity as x tends 
to 1 since the critical points of V (namely 1 and 02) merge to a single 
degenerate critical point. 
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4.3.1 The period function 

As in the case |a| < we can prove monotonicity of the period 

function T{a,l3,E) with respect to the energy E of the solution. The 
argument translates without much modification except for the fact that 
we want to prove that the Chicone test function N of (4.12) is positive on 
the interval [l,oo). 

Similarly, the period T{a, /3, E) goes to infinity as E ^ Hi. An exam¬ 
ple plot of T is given in Figure 4.5 (here H 2 — —0.4735 and Hi = —0.375). 



Figure 4.5: The minimal period as a function of the energy when a = 2 and 

(3 = 1. 


5 Counting solutions on X 52 


We are ready now to count the number of solutions of the Lichnerowicz 
equation on x 5^, bounded from above and away from zero, with g 
given by (2.2) and K of the form (2.3). 

If 0^/3^ > 4:R^(3^ there are no solutions. 

Otherwise there is always at least one constant solution, and we can 
rescale the metric so that R = + (3^. This scaling will be used in the 

remainder of this section. 

We have seen that when a = f3(V2 a^(3‘^/R^ = 4/3^, the only 

solution is (^ = 1. 

Suppose, next, that a < (3(^/2. We have seen that for 


T <Ti 


27r 


2V2n 
a// 32 - 2a2 


(5.1) 


the only solutions are constants (/> = 1 and (( = ((2 ■ 

Let Hi and H 2 be as defined at the beginning of Section 4.2. Now, for 
any (a, (3, E), with E S [Hi,H 2 ) the period function is continuous, strictly 
increasing in E, and satisfies 


T > Ti = lim T(a, (3, E). 

E\Hi 

The orbit with E = H 2 has infinite period, which implies that T{a,/3,E) 
tends to infinity as E tends to H 2 . It follows that for every T >Ti there 
exists precisely one value of E so that all solutions with energy E have 
minimal period T. Keeping in mind that the energy of the orbit is uniquely 
determined by the maximum value of (( at that orbit, for each value of E 
we obtain a one-parameter family of solutions, differing from each other 


12 










by the position of the maximum on the circle. From a geometric point 
of view these solutions can be considered to be identical, differing from 
each other by a translation along S^, which are isometries of g preserving 
the seed TT-tensor L. Here we will count the solutions modulo isometry, 
hence one solution for every energy level. 

As such, a solution on K with minimal period T = T{a, /3, E) provides a 
solution on M with period nT for any n S N* . Each such solution descends 
to X S'^ equipped with the metric gnT{a,/ 3 ,E) ,R=a^+p^- Set 




Figure 5.1: Plots of T, 2T, etc. as functions of energy at fixed a and j3 (first 
plot; compare Figures 4.4 and 4.5 for exact plots for specific values of a and 
/3), and as functions of (/>(0) (second plot). Rotating the plot by 7r/2 clockwise 
allows one to read the number of solutions as a function of T, as well as the 
corresponding values of </>(0), see Figure 6.1 below. 


To{a,(3) 


Ti{a,/3) ii\a\<(3/V2, 
T2{a,/3) ii\a\>l3/V2, 


(5.2) 


where Ti (resp. T 2 ) is given in Equation (4.11) (resp. (4.15)). Keeping 
in mind the constant solutions, our results so far can be summarised as 
follows (compare Figure 5.1): 


Theorem 5.1 Let a,/3>0, a ^ n G N. For any metric gT,a‘^+i}‘^ 

with T G {nTo{a, /3), (n + l)ro(a, /3)] there exist exactly n + 2 solutions of 
the Lichnerowicz equation modulo isometry. Two solutions are constant, 
and the remaining n solutions are not constant, and are invariant under 
rotations of the -factor. 


6 A bifurcation analysis 

One of the tools for constructing solutions of elliptic PDFs proceeds through 
bifurcation theory. The underlying idea is to study how the set of solu¬ 
tions evolves when varying a parameter. In this section we reexamine our 
problem from that point of view. 

We refer the reader to [7,23,25-27], [20, Sections 3.2 and 3.3], and 
references therein for a detailed account of this theory. For the reader’s 
convenience, we state the results that we will need in the paper. 

Consider a continuous mapping F : I x LI ^ B, where / is a non¬ 
empty interval of M, H is a subset of a Banach space A and B is another 
Banach space. We want to study the zero-level set of F, i.e. the set of 
pairs {X, x) G I X A for which F{X, x) = 0. 
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A pair (Aq, Xq) € / x 17 is called a bifurcation point for F if there exists 
a sequence (Ai)i>o converging to Aq with Xi € I and two sequences (a;^)i>o, 
(a;f)i>o in 17 such that 

• 'ii,xj ^ xf, 

• x-,xf Xo, 

• F(Xi,x^) = F{Xi,x‘l) = 0 . 

If the mapping is C^, this imposes that the differential dxF is not 
invertible at (Ao,a;o) since, otherwise, this would contradict the conclu¬ 
sion of the implicit function theorem. Two types of bifurcation points 
will be important for us in what follows. They are described in the next 
propositions. 

Proposition 6.1 (Fold bifurcation) Let (Ao,a;o) he a point in the zero set 
of F. Assume that dxF{Xo, xq) is Fredholm with index 0 and that the kernel 
of dxF{XQ,xo) has dimension 1 and is generated by Xi. Assume further 
that d\F(Xn, Xg) is not in the range of dxF{Xg,xg). Then there exists a 
neighborhood J x Lt' C I x Lt of (Ao,xo) and a C^-curve 7 :[/—)■ J x 17', 
where [/ C M is a neighborhood of 0, such that 

• 7(0) = (Ao,xo), 

• 7(0) = (0,a:i), 

• V(A, x) G J X 17', F{X, x) = 0 3t G U, (A, x) = 7 ( 7 ). 


For a proof of this proposition, we refer the reader to [23, Theorem 
2.3.1]. More information can be gained if we assume that F is C^: 

Proposition 6.2 (Fold bifurcation 2) Under the assumptions of Proposi¬ 
tion 6.1, there exists a linear form n G B*, n ^ 0, whose kernel is the range 
ofdxF{Xo,xo). Assuming further that F isC^ and pi{dxF{Xg,xg){xi,xi)) 
0 , then the curve 7 = ( 7 a, 7 x) C'" and 7 a( 0 ) = ■ 

In particular, upon shrinking the neighborhood J x LI' of (Ao,a;o), we 
have: 


• If > 0, then, for any X G J, 

#{xel7',F(A,a:)=0} 

• -(/7a(0) < 0, then, for any X G J, 

#{xe 17',F(A,x) =0} 


0 if X< Xq, 

1 if X = Xo, 

2 z/A>Ao. 


2 if X < Xo , 
1 if X = Xo, 
0 ifX>Xo. 


The second type of bifurcation was discovered in [7]: 

Proposition 6.3 (Pitchfork bifurcation) Assume that F is and that 
7 = (7 a, 7 a;) : U ^ I X LI is a curve of solutions: 

yt G U, F(7A(7),7a:(7)) = 0, 


with U a neighborhood of Aq in M such that 7a(7) = t, jxiXo) = xq- 
Assume further that dxF(Xo, xq) has a 1-dimensional kernel spanned by 
v G A and that D‘^F{Xo,xo){j{0),v) 7 ^ 0. Then (Ao,xo) is a bifurcation 
point for F and there exists a neighborhood J x LI' of (Ao,Xo) such that 
the set of solutions of F{X, x) = 0 consists of the union of two curves 
(7 and another one) intersecting (transversally) only at (Ao,a:o)- 
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6.1 T as a bifurcation parameter 

In this section we sketch the analysis of a bifurcation problem where T is 
considered as a bifurcation parameter. A similar more detailed presenta¬ 
tion of a bifurcation analysis, where a is the bifurcation parameter, will 
be given in Section 6.2. 

Consider a stable constant solution of (4.2), then candidate bifurcate 
solutions appear when the linearization of (4.2), namely 


(27r)^ d^v 

X ’2 


(PV 


{(pc)v, 


( 6 . 1 ) 


with suitable boundary conditions, has non-trivial kernel. This will be the 
case if and only if 


T 

2tt 


dW 

d(p 


(</>c) e N*. 


Indeed, note that solutions of (6.1), when they exist, come in two-dimensional 
families, parameterised e.g. by u(0) and dv /To set up a bifurcation 
theory argument with one dimensional kernel, one can consider those ele¬ 
ments of the kernel for which either 'c'(O) = 0, or u(0) = 0. One then finds 
that at each such value of parameters a new branch of solutions appears. 
This leads again to a picture as in the right Figure 5.1, at least near the 
intersection of the bifurcating solutions with the axis of constant solutions, 
except that now one finds apparently twice as many solutions. The reso- 



Figure 6.1; A bifurcation diagram for solutions, with 0(0) plotted as a function 
of T; compare Figure 5.1. 

lution of this apparent paradox is that two solutions with 0(0) = 1 which 
differ by the sign of d^(j)(fS) correspond to different solutions in the bifur¬ 
cation picture, while they were identified in our previous analysis: indeed, 
they have the same energy, and one can be obtained from the other by the 
isometry 0 i—>■ — 0. 

As should be clear from our previous analysis, the solutions on different 
bifurcation branches are actually the same solutions when the argument is 
allowed to run over K, but are interpreted as having a different periodicity. 

6.2 a as a bifurcation parameter 

In this section we study the equation (2.4) from the point of view of bifur¬ 
cation theory fixing g (i.e. T and i?), /3 and making a vary. In particular. 
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we no longer impose the conformal gauge R = From Theorem 

3.1, (j) depends at most on 'll;, so (2.4) reduces to 


As seen in Proposition 2.1, there is no solution to (6.2) if 

and a unique solution when ■= 2 ^R^ which is constant: 


4> = 4>o = 



In accordance with the terminology of [23, Remark 2.3.2], the point ip = (j>o 
is a subcritical fold bifurcation. For lower values of a, we get two branches 
of constant solutions going down to a = 0: 


p = p+ia) 

< 


R 

I 

(Nia) 

21/3^2 \ 

3/32 

3/32 

21/3 

N{a) J 

R 

1 

/,21/3R2 

1 

3/32 

3/32 

V /V(a) 

3 21/3 



where j = (—1 + i\/3)j2 and N is given by 


N{a) 



27a2/j4 _p 373 ^- 4 ^ 30 , 2^4 27a^P^ 


1/3 


Note that is a complex number. It can be shown that are real with 
0 < ^_ < (/)+. A plot of these solutions is given in Figure 6.2. From 



Figure 6.2: A plot of 4>+{a) (blue) and (p-ia) (red) with /3 = 1 and R = 
as a function of a. 


the shape of the potential, we see that p- is unstable while </>+ is stable 
within the class of solutions defined on intervals of K. To find potential 
bifurcation points, we look for p’s solving (6.2) such that the linearization 
of (6.2), namely 


(27r)^ 

T2 


(5/327 - - ^) 7 


(6.3) 


16 














admits a non-trivial solution We introduce the following function spaces: 

C'even('S'^;I^) := {? € C'^(S'^,M) is an even function of t/’} . 

These function spaces will be important to suppress the -translation- 
invariance of the set of solutions. We assume first that (j) is constant. 

Proposition 6.4 Bifurcations on the curve </> = ^+(0) occur for the fol¬ 
lowing values of a: 


a/c,± 



(6.4) 


where fc G N zs such that < "f • The values q;o,± = ±cimax 

correspond to fold bifurcations described earlier, while ak with k > 1 cor¬ 
respond to pitchfork bifurcations d la Crandall-Rabinowitz [7], There are 
no bifurcations on the curve (j) = (j)- (a). 


We note that the values ak,± can be rewritten as 

I ^ s 

2 


ttfc.i = ± 


3v^/32 


27r 


T 


i?3 - 2 ^ - 3i? ^ k ^ 


2tt 


T 


from which it follows that all values ak,± lie in the range [— Omax, cimax]- 


Proof: Since (f is constant, the right hand side of (6.3) is constant. Since 
f is 27r-periodic, this imposes the condition 

l(5/3V-7aV-'-B)=t"l^, (6.5) 

for some fc G N. The corresponding solution f is then, up to multiplication 
by a constant, 

C = cosikifj - ■0o)) . 

Values of a and (j) for which (6.2) and (6.5) hold can be found as follows: 
We introduce the polynomials 

P{X) = - i + a^ - RX^'^ , 

Q{X) = k^X"^ + i (^RX^ + 7a^ - 513‘^X^^ , 

which are obtained, for P, by multiplying the right hand side of (6.2) by 
(jf^ and setting X = and similarly for Q by multiplying (6.5) by (j>^ and 
setting X = The resultant of P and Q is given by 


0^/32 

4096 


+ 12 k^R - 4i?3 + 27a2/?4 


It is zero when a = 0 or when a = ak,± (see (6.4)). This means that when 
a = ak,±, P and Q have a common root given by 
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( 6 . 6 ) 



This value of X corresponds to (j)k ■= xl^'^ = 
It can be checked that 


so, for all values of A: > 0, (f>k is a stable local minimum for V. This 
proves that the bifurcation points along both branches (j)±{a) of constant 
solutions are located only on the curve (j) = (j)+{a). 

We now check that [7, Theorem 1] applies in this case. As we did in 
Section 4, to get rid of the S'^-invariance, we restrict the space of solutions 
to the Banach space C'|ygn(S'^,K) and restrict ourselves to the study of 
solutions (p to (6.2) belonging to this space. This restriction is actually not 
important since any solution p to (6.2) admits a point i/'o where p'{ipo) = 0. 
It follows from the Cauchy-Lipschitz theorem that (j){'ijjo+Sip) = 
yS'tp € M. Translating the solution, we can assume that p G C'|ven('S'^) I^)- 

We let 


F : C2ven(5^R) n {((. > 0} X M ^ 


be the following operator: 


F{p,a) 


(27r)^ (Pp 

T 2 


^ Fa^P - Rp^ ■ 


At points {pk, o:k,±), the linearization of F has a 2-dimensional kernel 
generated by the following two vectors 

vi := i6pi,Sai) = {-2ak,±pk, - Rpl - 7(^1 ,±). 

V 2 '■= iSp 2 , Sa 2 ) = (cos(A:'i/')) 0) • 


The derivative 

has one-dimensional kernel generated by dp 2 = cos{k^p), and its image is 
the kernel of the map 


f i~> / f{'p)cos{k'p)dp. 

Js^ 

This is the reason why we restrict to the space of even functions, otherwise 
the kernel of D^F(pk,ak,±) would be two-dimensional, similarly for the 
cokernel, thus failing to satisfy the assumptions of [7, Theorem I]. The 
only condition that remains to be verified is that 

F''{pk,ak,±){vi,V2) ^ R{F'{pk,ak,±)). 


This actually follows from a straightforward calculation: 


F"{pk,ak,±){vi,V2) = 


ak,± 

4 


^ - ^bptP"^ ) cos{kp) 
Vk 

2 


142 


27r 

f 


kr + \2\R cos{kp). 
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□ 

Our next step is to obtain a better understanding of the curves of non¬ 
constant solutions. To label the branches solutions, we define the index 
of a solution. Given a non-constant solution (j), we have, for all -ip G S^, 
, (p{ip)) ^ ((/)+(«), 0). So a non-constant solution ^ is a curve in 
\ {(</)+(«), 0)}. We define its index as the class of </> in 

This index is constant along a curve of solutions, except at the bifurca¬ 
tion points on the curve a i—)■ (p-^. (a) where the index is not defined. Each 
solution on the curve a i—>■ 4>-{a) has index zero while each bifurcation 
point (/)+(Q;fc^±)) is the limit point of two curves of non-constant 

solutions with index k: 

Proposition 6.5 For all k > 1 such that § there exist two 

curves 

cifc,+) 4^k,±io:) € C'|ygn(5'^, M) 

of solutions to (6.2) of index k which are 2'!r/k-periodic. The curves are 
obtained one from the other as follows: 

(pk-W = (l)k,-\- (V'+ • 

No bifurcations occur on these curves except at the points ak,±. These 
solutions, together with the solutions lying on the curves a i—>■ (j)±{a), 
exhaust the set of solutions to (6.2). 


Before proving this proposition, we need the following lemma: 

Lemma 6.6 The period T{a, E) of the solutions to (6.2) with energy E de¬ 
pends analytically on{a,E)foraG (—dmax, cimax) and E G iy{4>+{a)) ,V{(f>-{a))). 


Proof: The proof is based on a rewriting of (4.9): 


T{a,E) = V2 


^dmax(a,iJ) 
' <^min(a.E) 


df) 

^E-V{cP,a) ’ 


where < 4>max(pi,E) are the two solutions to V{(f>,a) = if in 

the range {(j)-{a) ,oo). Note that since 


d(j) 


(0min(o:5 E'j , if) , 




i{a,E) ,E)^0, 


the analytic implicit function theorem shows that (prnin and ^max are an¬ 
alytic functions in a and if. Given and Eq satisfying the assumptions of 
the lemma, we choose an arbitrary value (po G {(pminicto, Eq) , ^max(Q;o, Eq)). 
Given {a,E) close to {ao,EQ), we split (4.9) as follows: 


T(a,E) = 


rf’o 


d(p 




d(p 


U^Ua,E) y/E -V{(p, a) 


'do 


y/E-V{eP,a) 
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We show how to rewrite the first integral so that its analyticity in the 
vicinity of (ao,-Eo) becomes apparent. Note that 


E -V{(j),a) = V ((/)inin(a, E),a)-V {(j>, a) 

dv 


f" dV 

= -{(j> - (j)n,in{a,E)) / —{X(j>+{l-X)(j)^inia,E))dX 
Jo 

=i 4>ni\n{.a,E) ,a)dX, 

Jo 0(j) 


where we set 4> = xj + ^min(a, E). So, 


c4>a 


d(j) 


4„i„(a.£;) yjE -V{(j),a) 

ry/ 4>o—4>niiii{<^,E) 


= 2 


dx- 


\j~ fo + «^min(a, E) , a)dX 

A similar rewriting of the second integral yields 

d(j) 


4o \/E-V{(j),a) 

f 

= 2 


Cy/ dmax(<^,£') 00 


dx-i 


>0 


\//o^ + </’max(a, E) , a)dX 


where (j) = (/)max(Q!, E) — x\. The function 

{E,a,X-)^-f ^iXx'^_ + 4>minia,E) ,a)dX 
Jo d(j) 

is clearly analytic and positive for all x+ € 0, y/</>ma.x(o:, E) — 0o 

dV 


since 


d(j) 


(^min(^j E^ , o) < 0 . 


This is enough to conclude that 

f*° d(j) 

J4>^Ua,E) y/E -V{4),a) 

ry/ 4>o~4>niiii{<y,E) 
= 2 / 


fix_ 


fo + ^min(a, E) , a)dX 

is analytic in (a, if) in a neighborhood of (ooi^o)- Similar arguments 
apply for the second integral. □ 


Proof of Proposition 6.5. We first remark that since the energy iJ 
defined in (4.3) is conserved, solutions ^(V') to (6.2) with index k are 
actually periodic with minimal period 27r/fc. We select a non-constant 
solution {ao,4>o) with index k and energy E^g- Since the index is locally 
constant, all solutions nearby, potentially with a different a, are 270Jk- 
periodic. 
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From the analysis in Section 4, for all values of a between zbamax 
the period ^ r{E) of a solution (/> with energy E = (p) is strictly 

increasing with respect to E. Since we are restricting ourselves to solutions 
belonging to C'|ven('S'^)®): we have ^(0) = 0 so i? = fo((/)(0)). Note that 
since the derivative of T with respect to E is strictly positive, for all a 
near ag there exists a unique value E^ € (fo(0_i_(a)), fo((/)_(a))) of the 
energy so that the solution with energy E^ has period 27r/fc. E^ depends 
smoothly on a by Lemma 6.6. We let denote the unique solution 

(f) > (/'+(q;) to V{<j)) = V{<j)-{a)). From the shape of the potential V, there 
exist exactly two values ipminict, k), </>max(a, k) so that 

(l)-{a) < k) < ^+(q;) < (/'max(Q!, k) < (j)*_{a), 


and 

^ (0min(c^5 ^)) — fo(^max(Q^j ^)) — ■ 

These two values map smoothly to two solutions of (6.2). Thus we have 
proven that near a value ao for which there exists a solution (j)Q with index 
k, there exist two and only two distinct curves of solutions with index k. 

Note that a 27r/fc-periodic solution goes from k) to (/)max(cr, A:) 

in an interval of length ir/k and then goes down from ^max(Q;, k) to ipminict, k) 
in the same amount of time. Hence, translating the solution (j) with 
(()(0) = 4>minio:,k) by ir/k we get the solution (j){0) = 4>max{oi,k) and vice 
versa. 

Let 

4 c M 

denote the set of values for which there exists a pair of solutions of index 
k. The previous analysis shows that 4 is an open subset. Assume that 
4 contains a boundary point a^o which is not ak,±, Let Ui G 4 be such 
that ai -G ttoo- The corresponding functions 4>i with period 27r/fc all have 
€ {(j)-{ai) ,(j)*_{ai)). Without loss of generality, we can assume that 
(t>i{0) converges to some limit ^oo(O) € [(j)-{aoo) ,^l(aoo)]- ^^'oo(O) cannot 
be (j)-(oioo) nor 4>*_{aoo) since the period of the functions 4>iS would grow 
unbounded. 

If (/)oo(0) 7 ^ ^+(aoo), by continuity of the period with respect to initial 
data, the solution /poo to (6.2) is periodic with period 2tt jk. The previ¬ 
ous argument shows that aoo is an interior point of 4, a contradiction. 
Thus the only endpoints of Ik are on the curve p = </>+(«), this is to say 
bifurcation points we found in Proposition 6.4. 

The question now arises, whether new solutions occur with values of a 
larger than a^, or smaller, or both. We will see that, for all k such that 

4 contains an interval of the form {ak,+ — e,ak,+)- Since 

(6.2) only depends on Ik also contains the interval {ak,-,ak,- + e)- 
We let 

(- 5 ,( 5 ) 9 1 1-9 {a(t), pkit)) 

denote a differentiable curve of non-constant solutions passing through 
{ak,+ i4‘k) at the value t = 0 of parameter t, where 
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compare (6.6). The existence of this curve has been established in Propo¬ 
sition 6.4. We expand a and (j)k in terms of the parameter t as follows: 


a{t) = ak,+ -I- , 

4>k{t) = 4>k + t$l. + (pk + t^Pk + o{t^), 


where (l>\ = cos^kip) and insert this development in (6.2). From the terms 
linear in t in (6.2), it follows that dj, _|_ =0. Looking at terms of order 
we find that = Acos(2fc^/>) -|- ecos{kip) -|- where 


A = 

12fc2 



and 


61/V/5 



There is no loss of generality in assuming that e = 0 since this can be 
reabsorbed in the definition of t. More importantly, /i depends on oi^ 
this is why we need to consider terms cubic in t in (6.2). The expression 
for is obtained by setting to zero the coefficient of cos{ktp) in (6.2): 



The numerator of this expression is decreasing on [0,oo) when seen as a 
functmn of k so it is bounded from below from the value it takes when 
which is 5R/2. Since we have a = ak,+ + -I- 0{t^) 

with _(_ < 0, we deduce that a{t) < ak,+ for small values of the param¬ 
eter t. This concludes the proof of the claim. 

Ik being connected with endpoints afc,±, we conclude that Ik = {ak,-,ak,+)- 

□ 


An illustration of the last two propositions is given in Figure 6.3. 
Summarising, we have proved: 

Theorem 6.7 Let T,R G M"*", and let (3 > 0 be defined in (2.5). Consider 
a metrie g = gf and define /Cmax to he the largest integer k such that 


2 



< R. 


Depending on the value o/a G M, Equation (6.2) has, up to translation in 
the - direction, 
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Figure 6.3: An illustration of Propositions 6.4 and 6.5 with T = 27r, R = 33 
and /3 = 3. On all plots a varies along the horizontal axis. On the first line, 
the plot on the left shows the three curves of (magenta), (/>+(a) (marine 

blue) and (p*_{a) (yellow). Dots indicate the position of the bifurcation points. 
Where curves merge, there are actually two points which almost coincide. One 
corresponds to the fold bifurcation while the second one is a pitchfork bifur¬ 
cation. All other points are pitchfork bifurcations. The second plot shows the 
value at the origin of the solutions of index 1 (brown), 2 (red), 3 (green) and 4 
(blue). And the third plot is a zoom of the second one near amax- The graphs 
on the second line show the energy H(cj), (p) of the solutions. 


• no solutions if a ^ [-Omax, Omax], where a^ax = 

• only one solution if a = ittmax) 

• two constant solutions and k non constant SO (3)-symmetric solu¬ 
tions of index 1, k when a G {ak,-,cxk+i,-\ U [afe+i_+,afe,+) if 
k < fcmax or when a G {ak-,ak,+) if k = /c^ax- 

For further reference we note the following 

Proposition 6.8 Let denote the linearization of the Lichnerowicz equa¬ 
tion at 4>. Then 

1. If (j) is a constant solution then has no kernel except at the bi¬ 
furcation points described above, where the dimension of the kernel 
equals two. 

2. If (j) is one of the non-constant solutions above, then has a non¬ 
trivial kernel. 


Proof: 1. Let 

L^v={2^g + V''{(P))v = G, ( 6 . 8 ) 

where (/) is a constant solution. Let (pt be an eigenfunction of the Laplace 
operator on S 2 with eigenvalue —t(l+ 1), £ G N. Set vi = fpi,v)si, where 
(•, •)52 is the standard L^-product on S'^. Equation (6.8) implies 


(27r)^ dfve 
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(6.9) 


Thus vi = Acos{k'tp) + i?sin(fc^), with 

^^e = {v"{<p) + i{e + i)). 

1 ^ 

This equation together with the condition V'{(j)) = 0 gives two polyno¬ 
mial equations for (j). A calculation shows that the resultant of these two 
equations has no real roots ior £ > 1, which establishes the result. 

2. Let (() be a non-constant solution of the Lichnerowicz equation, then 
d(j>/d'ip is a non-trivial element of the kernel. □ 


7 CMC slicings of Nariai and Schwarzschild 
- de Sitter 


The initial data sets constructed above are invariant under rotations of S'^. 

It follows from the generalised Birkhoff theorem [8, 28, 30] that the associ¬ 
ated maximal globally hyperbolic developments are subsets of Schwarzschild- 
de Sitter space-time or Nariai space-time. (Here the de Sitter solution is 
considered as being a member of the Schwarzschild-de Sitter family.) As 
such, we have thus been constructing the geometry of CMC slices in the 
Schwarzschild-de Sitter and Nariai space-times; compare [2,10], where the 
constants K and C are given in our notation by 


K = T , 


C = 


2a 

\/3A3/2 ■ 


7.1 Nariai 


We start by recalling the standard form of the Nariai metrics, 

1 

g = _(A-Ar2)ci^2 + ___ + _d^7^ (7.1) 


with A € M. In fact, rescalings of the t- and r-coordinates allow us to 
achieve A € {0, ±1}. 

For the purpose of the analysis that follows, the key property is that 
the angular sector of the metric (7.1) is both t- and r-independent. 

Suppose that our initial data set (M = x S‘^,g,K) arises from a 
periodic hypersurface in Nariai space-time. This is compatible with (2.2) 
if and only if 



(7.2) 


In particular (j) must be constant. 

Next, the extrinsic curvature of any spherically symmetric hypersur¬ 
face in a Nariai space-time will have trivial components in the spherical 
directions. It follows that (2.3) is compatible with the Nariai form of the 
metric if and only if 


^ 3 ^ ^ J 


(7.3) 


Equations (7.2)-(7.3) give necessary conditions for existence of an em¬ 
bedding of our initial data sets into a Nariai space-time. 
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Let us show that these conditions are sufficient. For this, it turns 
out that we only need to consider the region where r is a time-function. 
Obvious renamings bring (7.1) to the form 

1 

g = ^ ~ X)dx^ + — A > 0. (7.4) 


Let n = \/Afi — Xdt denote the held of unit normals to the slices t = const. 
The mean extrinsic curvature T^it) of those slices is readily calculated to 
be 

TN { t ) = = 


\/det^ 


-dai^Jdei= dt{V- A) 


At 


If (7.2) and (7.3) hold, we can hnd an embedding of our initial data set 
into a Nariai space-time by hnding values of A and t so that TAr(t) = r: 

At 


A{A-T^)t^ = -Ar^ 


(7.6) 


VAt2 - A 

Hence, the development of initial data satisfying (7.2) and (7.3) corre¬ 
sponds to a Nariai space-time with 

1 if t2 > A, 

0 if r2 = A, 

-1 ifT2<A. 


A = 


It is instructive to locate the solution corresponding to the Nariai space- 
time on the graphs of Figure 6.2. A good indication of the location of this 
solution is given by the stability of the point it corresponds to. Without 
assuming the normalization R = , the potential V associated to 

the ODE (6.2) is given by 


VW 



Using the value of j3 from (2.5) and the values of a and (j) given by (7.2) 
and (7.3), we hnd 

V"W = ^{A-r^). 

This means that, if a constant solution to the constraint equations cor¬ 
responds to the Nariai space-time with A = —1 (resp. A = 0, A = -fl), 
the corresponding point on Figure 6.2 is on the upper branch (resp. the 
rightmost point of the diagram, the lower branch) because the value of (j) 
corresponds to a stable (resp. degenerate, instable) critical point of V. 
We can rewrite the value of a given by (7.3) as follows: 


a = 



3,12 ■ 


Letting vary in the range [0, oo) at f3 and R hxed, we see that a increases 
from 0 for = 0 to Omax when ^ and then decreases to 0 when 

> oo. We conclude that, hxing (3 and R and allowing a to vary, the 
“Nariai point” (a, </>) can lie anywhere on the constant solution curves, 
such as the curve in Figure 6.2, their location being determined by the 
relative value of r and A. 
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7.2 Schwarzschild-de Sitter 


Whenever (7.2) and (7.3) do not hold, in particular when (p is not constant, 
the development of the initial data will be a Schwarzschild-de Sitter space- 
time. 

We shall (essentially) use the Hawking mass to determine the corre¬ 
sponding mass parameter. 

For this, we continue with the following calculation, somewhat more 
general than needed: Consider an (n -I- l)-dimensional metric, n > 3, of 
the form 

g = —fdt^ -\—— + hAB{x'")dx^dx^ , (7.7) 

/ ''-V-^ 

— :h 

where h is a Riemannian metric on a compact manifold M with constant 
scalar curvature R] we denote by x^ local coordinates on M. As discussed 
in [3], for any to G K and £ G K* the function 


/ 


R 

(n — l)(n — 2) 


2 to 

j.n—2 £2 


(7.8) 


leads to a vacuum metric. 


R^iu — ^2 9^^v 1 

thus f is a constant related to the cosmological constant as 

1 _ 2A 
P n{n — 1) 


(7.9) 


(7.10) 


We note that there are five different families of Schwarzschild — deSitter 
space-times, with distinct global structure, depending upon the values of 
i? > 0, TO G M and A > 0: 

1. If TO < 0, / is strictly decreasing from oo to —oo on the interval 
(0, oo). Thus there exists a unique value Vc such that f{rc) = 0. The 
region r < Xc describes a static space-time with a naked singularity. 
We call these space-times negative mass Schwarzschild - de Sitter. 

2. If TO = 0 we are in de Sitter space-time, with a single first-order zero 
of /. 

3. If 0 < 6-\/2 Ato < the lapse function / vanishes for two distinct 
positive values of r which we denote by r_ and r_|_. This corresponds 
to the “usual” Schwarzschild - de Sitter space-time which we call 
subcritical Schwarzschild - de Sitter. 

4. If the lapse function / has a single double zero, we 

refer to this situation as extreme Schwarzschild - de Sitter. 

5. For larger values of to, / remains negative on the whole interval 
(0,oo). The singularity at r = 0 corresponds to a big bang-type 
singularity, without any horizons. These space-times will be referred 
to as supercritical Schwarzschild - de Sitter. 

The value of to which separates cases 3 and 5 can be found as follows. The 
polynomial rf, where / is given by (7.8), is a third order polynomial: 


i? Ar3 

rt = —r — 2m -—. 

2 3 
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When m > 0, rf always has a unique negative root since / increases 
strictly on the interval (— 00 ,0) from —00 to 00 . The number of real roots 
of rf is given by the sign of its discriminant: 

A = ^(^R^- 72m^A) . 

When A > 0, r/ has three real roots, so, from the discussion above, two 
positive roots. This corresponds to the subcritical case. When A < 0, 
rf has only one real root which has to be the negative one. So / keeps 
constant negative sign on (0,oo). This corresponds to the supercritical 
case. 

To determine the space-time associated with our initial data set we 
consider the surfaces Ma,b '■= {r = a,t = b} C M and associate to them 
the following integral: 




[ 0+0-, 
\Mr,t\ JMr.t 


where is the area of with 


0+ =g^^KAB±H, 


(7.11) 


(7.12) 


where H is the mean extrinsic curvature of Mr.t within the slice of constant 
time. And, as usual, Kij is the extrinsic curvature tensor of the slices of 
constant time. The integral is closely related to the Hawking mass of the 
manifolds Ma,6- 

The key fact is that, while each of 0+ and 0- depends on choices made, 
their product does not, and therefore I{Mr,t) is an invariant determined 
solely by Mr,ti which justifies the notation. Note that in our case is 
proportional to the area of Mr,t, in fact 

r"^ = 4>‘^. (7.13) 


We will work in the region where / > 0. Using the time function t 
of (7.7) we have Kij = 0. The metric induced on the slices of constant 
time is f-^dr'^ + r'^h. The field, say b, of unit normals to the level sets of 
r reads v = ^/fdr- Denoting by D the covariant-derivative of the metric 
gijdx^dx^ induced on the level sets of t, we hnd 


H 


hence 


Dkiy’' = dki^/detg^jiy'') 

^det g^j 

(n- l)^/f 

r 




(7.14) 




(n-l)Vr-" 
A (R 2m 


Remaining in dimension n = 3, we return to our initial data set 
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fi,K) with K given by (2.3), thus 


9 = 

K = 



(7.16) 

(7.17) 


The field of unit normals to the level sets of i/' is = '^<j> ^9^, leading to 
a mean curvature 

H = 

Using 


1 _ 

4(n — l)7r9^(() 

f(j)^ 


-6 


ABj^ o I ‘ 


we conclude that, in dimension n = 3 

2 


I{Mt,r) = 


( 3TTd^4>\ 


V ) 

From (4.3) at energy E we have 


1 ( 2tt d<f> 


R 


_4( 

3 V6 


= E 


T did J 16 


Inserting into (7.20) one finds 


48 


48 


(7.18) 


(7.19) 


(7.20) 


(7.21) 




a(j) 


-6 




-4 — - 
9 


T 

3 " “^/6 

^ 2Ta4>~^ 


3V6 


= n&E 


Ara \ 2R 2/3^ 4t^ 

+ ^ 9“ 


(7.22) 


Comparing with (7.15), we conclude that we have constructed initial data 
on a CMC slice in a Schwarzschild-de Sitter space-time with 


m 4U . (7.23) 

As such, this formula still holds in the regions where / < 0 when 
appropriately understood. Indeed, note the following ambiguity: In the 
region where / > 0, the expansions have been calculated with respect to 
the vectors 


e-h = -^dt + \/7'9r , e_ = -^dt - y/fdr 


These vectors satisfy g{dt,e±) = < 0, which implies that they 

have the same time orientation. If one wishes to preserve this property in 
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the region where / < 0, where dr is time-like, one should use instead e.g. 
the following pair of lightlike vectors: 


e-h = V\T\dr - ’ e- = y/\T\dT + ' 

It is natural to inspect Figure 6.3 from the point of view of the CMC slices 
in the Schwarzschild-de Sitter space. In Figure 7.1 we show representative 
plots of the mass of the solutions. The value of the mass corresponding to 



Figure 7.1; Masses of the Schwarzschild-de Sitter space-times as functions of 

а. The values of T = 27r, R = 33 and /3 = 3 coincide with the ones in Figure 

б. 3. We chose r = 5 and A = 77/6, consistently with (2.5). The plot on the left 
shows the masses associated to the constant solutions to (6.2). The middle plot 
includes further all non-constant solutions to (6.2). The third plot is a zoom in 
the region a ~ Omax- The cross symbol indicates the values of the parameters 
corresponding to the Nariai space-time. 


the cusp on the right (i.e. when a = Omax) can be computed explicitly in 
terms of P and r: 


m = 


27^2 


( 2 ^ 3 ^ 


— r 


Similarly for the left cusp: 


y2.R3/2 

27^2 


( 2 ^ 3 ^ 


-I- T 


Hence, choosing r such that > 12/32 we obtain solutions with negative 
masses, see Figure 7.2. 



Figure 7.2: An example of a family of solutions to the Lichnerowicz equation 
containing Schwarzschild-de Sitter space-times with negative mass. Here, T = 
27r, R = 33, /3 = 1 and r = 4 (hence, A = 35/6). The plot shows the masses of 
the constant solutions to (6.2) as a function of a. 
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So far we have identified which Schwarzschild-de Sitter space-time will 
arise from our data. It appears of interest to enquire where the initial data 
set will lie in the associated space-time. We note the following: 

Consider, first, the data set determined by constant-^ solutions, (p = 
In the case that the time development is a Schwarzschild-de Sitter 
space with positive subcritical mass, let 0 < r_ < r+ < oo denote the area 
radius of the horizons in the space-time with the mass determined from 
the parameters of the solution, as described above. 

From (7.13), these constant solutions correspond to constant r hyper¬ 
surfaces. Since the induced metric (pP^g is spacelike, it should be clear from 
the geometry of the Schwarzschild - de Sitter space-time that either the 
initial data set embeds in a Nariai space-time, or the CMC surface lies 
inside the “hole horizon”, whether white or black, that is in the region 
{r < r_}, or above the cosmological horizon, in the region {r > r+j. 
In the Schwarzschild - de Sitter case solutions with nearby energy will 
remain in the same region, without crossing any horizons. As the en¬ 
ergy is increased, the solution oscillates between a minimum smaller than 
4>+{a) and a maximum larger than (j)^(a). This means that the solution 
keeps intersecting the original region determined by </>+(«). Now, a space¬ 
like hypersurface cannot cross a connected component of the collection of 
horizons back and forth. Hence the whole hypersurface must be entirely 
contained in the original region. We conclude that, for Schwarzschild - de 
Sitter space-times with fixed positive subcritical mass, 

1. All initial data sets which can be continuously deformed, by changing 
the parameters or the energy, to a constant solution lying under a 
hole horizon are entirely contained under that hole horizon, and 

2. all initial data sets which can be continuously deformed, by changing 
the parameters or the energy, to a constant solution lying above the 
cosmological horizon are entirely contained above the cosmological 
horizon. 

We conclude that the collection of complete periodic CMC hypersur¬ 
faces in Schwarzschild - de Sitter with fixed positive subcritical mass has 
at least three distinct components. Allowing the mass to vary, solutions 
belonging to the boundaries of the components have either zero or extreme 
mass. The latter possibility is illustrated in Figures 7.3 and 7.4. 

An interesting fact one might guess from Figures 7.3 and 7.4 is that 
there does not exist any non-constant solution lying entirely inside the 
hole region. This fact is proven in [10, Section 3.2]. 
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